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Abstract. We present an exact solution to the A"~ -proton bound state problem formulated in the mo- 
mentum space. The Is level characteristics of the kaonic hydrogen are computed simultaneously with the 
available low energy K~p data. In the strong interaction sector the meson-baryon interactions are de- 
scribed by means of an effective (chirally motivated) separable potential and its parameters are fitted to 
the experimental data. 

PACS. 11.80.Gw Multichannel scattering - 12.39.Fe Chiral Lagrangians - 13.75.Jz Kaon-baryon interac- 
tions - 36.10.Gv Mesonic atoms and molecules, hyperonic atoms and molecules 
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1 Introduction 

We developed a precise method of computing the meson- 
nuclear bound states in momentum space. The method 
was already applied to pionic atoms and its multichannel 
version was used to calculate the Is level characteristics of 
pionic hydrogen , ]J . In the present work we aim at simul- 
taneous description of both the Is level kaonic bound state 
and the available experimental data for the K~p initiated 
processes. 

Until recently the old Deser-Trueman formula ^ was 
used to determine the strong interaction energy shift and 
width (of the Is level) in kaonic hydrogen from the K~p 
scattering length and vice versa. Recently, the Deser - 
Trueman relation was modified to include the isospin ef- 
fects and electromagnetic corrections [3l. Our exact so- 
lution of the K~p bound state problem allows to check 
the precision and limitations of those approximate ap- 
proaches. However, one should not forget that the strong 
interaction part of the scattering length is not a directly 
measured quantity and its determination from the scat- 
tering data is always model dependent. 

The treatment of the kaon-nucleon interaction at low 
energies requires a special care. Unlike the pion-nuclcon 
interaction the KN dynamics is strongly influenced by the 
existence of the yl(1405) resonance, just below the K'p 
threshold. This means that the standard chiral perturba- 
tion theory is not applicable in this region. Fortunately, 
one can use non-perturbative coupled channel techniques 
to deal with the problem and generate the yl(1405) reso- 
nance dynamically. Such approach has proven quite useful 
and several authors have already applied it to various low 
energy meson-baryon processes ^-[Sj. Whether the recent 
experimental results on kaonic hydrogen from the DEAR 
collaboration [9] are consistent with the older KEK results 



[lOj and whether they fit into the picture drawn by the 
chiral models represents a question which is addressed by 
the theory [8[|TT] as well as by the coming SIDDHARTA 
experiment. 



2 Formalism 

Our approach to solving the meson-nuclear bound state 
problem in the presence of multiple coupled channels was 
given in Ref. Tj. Here we just remark that the method is 
based on the construction of the Jost matrix and involves 
the solution of the Lippman-Schwinger equation for the 
transition amplitudes between various channels. Bound 
states in a specific channel then correspond to zeros of 
the determinant of the Jost matrix at (or close to) the 
positive part of the imaginary axis in the complex mo- 
mentum plane. The zeros are computed iteratively and if 
only the point-like Coulomb potential is considered in the 
K~p channel the method reproduces the well known Bohr 
energy of the Is level with a precision better than 0.1 eV. 

We follow the approach of Ref. [4] when constructing 
the strong interaction part of the potential matrix. In this 
model the yl(1405) resonance is generated dynamically 
by solving coupled Lippman-Schwinger equations with in- 
put effective (chirally motivated) potentials. The reader 
should note that our approach differs from the recently 
more popular on-shell scheme based on the Bethe-Salpeter 
equation, unitarity relation for the inverse of the T-matrix 
and on the dimensional regularization of the scalar loop in- 
tegral [T5]. Further, while the authors of Ref. [4J restricted 
themselves only to the first six meson-baryon channels 
that are open at the KN threshold we employ all ten 
coupled meson-baryon channels: K~p, K'^n, tt^vI, tt+Z"^, 
ttOZ"", TT-S+, r]A, f]Z:", K+E-, and K"S°. 
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The strong interaction potentials are constructed in 
such a way, that in the Born approximation they give the 
same (up to 0{q'^)) s-wave scattering lengths as are those 
derived from the underlying chiral lagrangian. Here we use 
them in the separable form 



V^J{k,k') -~ 



1 M, 
2Ei uji 

1 



1 M, 



2Ej ujj 



(1) 



in which the momenta k and k' refer to the meson-baryon 
cm. system in the i and j channels, respectively, and the 
kinematical factors ^/WjJ{2EjUJj) guarantee a proper rel- 
ativistic flux normalization with Ej , Mj and ujj denoting 
the meson energy and the baryon mass and energy in the 
cm. system of channel j. The off shell form factors gj{k) 
introduce the inverse range radii aj that characterize the 
radius of interactions in various channels. Finally, the pa- 
rameter / stands for the pseudoscalar meson decay con- 
stant in the chiral limit and the coupling matrix Cy is 
determined by chiral SU(3) symmetry and includes terms 
up to the second order in the meson cm. kinetic energies. 
For the first six channels the couplings Cij were listed in 
[4] and we intend to publish the remaining coefficients in 
a more elaborate paper [T3]. For illustration, we show just 
the couphng of the elastic K~p process. 



CK-p,K-p — —Ek — 



El - ml 
2Mn 



3 ' 



2Mn 



+ Amliho + &o) - El{2dD + 2d„ + di) . (2) 

Here mx and Ex denote the kaon mass and energy in the 
center-of-mass frame, Mq stands for the baryon mass in 
the chiral limit, and the parameters F, D, and 6's and d's 
represent coupling constants that appear in the underly- 
ing chiral lagrangian (see [1] for more details). The origin 
and relevance of the various terms present in Eq. ^ was 
discussed thoroughly in Ref. [7]- In general, the coefficients 
Cij include contributions from the meson-baryon contact 
interactions as well as the direct and crossed Born terms. 
However, in contrast to [7^ our model is based on the static 
(heavy) nuclcon approximation adopted by the authors of 
Ref. [1] in which the underlying lagrangian is expressed in 
a fixed reference frame. 

The potential of Eq. ([T]) is used not only when solving 
the bound state problem but we also implement it in the 
standard Lippman-Schwinger equation and compute the 
low energy KN cross sections and branching ratios from 
the resulting transition amplitudes. 



3 KN data fits 

The parameters of the chiral lagrangian which enter the 
coefficients C'ij and the inverse range radii aj determining 
the off-shell behavior of the potentials are to be fitted 
to the experimental data. Before performing the fits we 
reduce the number of the fitted parameters in the following 



way. First, the axial couplings D and F were already fixed 
in the analysis of semileptonic hyperon decays [14j , D = 
0.80, F = 0.46 {gA ^ F + D ^ 1.26). Then, we fix the 
couplings bu and bp to satisfy the approximate Gell-Mann 
formulas for the baryon mass splittings. 



Ms - Mn = 



.{ml -ml) 



16, 



Me-Ma = -^bD(m 



(3) 



which gives bo = 0.064 GeV-^ and hp = -0.209 Ge\-^ . 
Similarly, we determine the coupling 6o and the baryon 
chiral mass Mq from the relations for the pion-nucleon 
sigma term g^^n and the proton mass, 

CTttW = -2ml{2bo + bD + bp) , 

Mp = Mo - Aml{bo + bD- bp) - 2ml{bo + 2&f) .(4) 

Since the value of the pion-nucleon a-term is not well de- 
termined we enforce four different options, ct^at = 20 — 50 
MeV, which cover the interval of the values considered 
by various authors. Finally, we reduce the number of the 
inverse ranges aj to only five: axN, olt^a-, cut^e, ctrjA/s 
(for both the t]A and rjS channels), and axE- This leaves 
us with 11 free parameters: the five inverse ranges, the 
meson-baryon chiral coupling /, and five more low energy 
constants from the second order chiral lagrangian denoted 
by da, dp, do, di, and d2- 

The fitted low energy KN data include the three pre- 
cisely measured threshold branching ratios [15] 



2.36 ±0.04, 



Rn 



R„ 



a{K^p TT^U^) 
a{K^p charged particles) 
a{K^p — *■ all) 
a{K-p n^A) 
a{K^p all neutral states) 



= 0.664 ±0.011, 



0.189 ±0.015,(5) 



and X^p-initiated total cross sections. For the later ones 
we consider only the experimental data taken at the kaon 
laboratory momenta p^AS = 110 MeV (for the K~p, K^n, 
TT+E-, 7r"r+ final states) and at plab = 200 MeV (for 
the same four channels plus tt^A and tt'^S'^). Our results 
show that the inclusion of the cross section data taken 
at other kaon momenta is not necessary since the fit at 
just 1 — 2 points fixes the cross section magnitude and 
the energy dependence is reproduced nicely by the model. 
Finally, we include the DEAR results [9] on the strong 
interaction shift AE^ and the width F of the Is level in 
kaonic hydrogen: 

AEn{Is) = (193±43) eV, F{ls) = (249±150) eV . (6) 

Thus, we end up with a total of 15 data points in our fits. 

Our results are summarized in Tables [T][31 The first ta- 
ble shows the results of our fits compared with the rel- 
evant experimental data. The resulting per data point 
indicate satisfactory fits. It is worth noting that their qual- 
ity and the computed values do not depend much on the 
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Table 1. The fitted KN threshold data 



a^N [McV] 




AEn [eV] 


r [eV] 


7 


Rc 


Rn 


20 


1.33 


232 


725 


2.366 


0.657 


0.191 


30 


1.36 


262 


697 


2.365 


0.657 


0.190 


40 


1.37 


253 


710 


2.370 


0.657 


0.189 


50 


1.40 


266 


708 


2.370 


0.658 


0.190 


exp 




193(43) 


249(150) 


2.36(4) 


0.664(11) 


0.189(15) 



Table 2. Chiral lagrangian parameters (6o and d's in 1/GeV): 



a^N [MeV] 


bo 


Mo [MeV] 




/ [MeV] 


do 


do 


dp 


di 


d2 


20 


-0.190 


997 


-0.016 


108.6 


-0.385 


-0.368 


-0.817 


0.396 


0.152 


30 


-0.321 


864 


0.001 


100.0 


-0.354 


-0.206 


-0.522 


0.406 


-0.211 


40 


-0.453 


729 


0.006 


108.9 


-0.484 


-0.151 


-0.459 


0.448 


-0.280 


50 


-0.584 


594 


0.007 


108.8 


-0.747 


-0.092 


-0.429 


0.567 


-0.349 


27 [4] 


-0.279 


910 


-0.002 


94.5 


-0.40 


-0.24 


-0.43 


0.28 


-0.62 



exact value of the a.,rN term. Tables [2] and [3] show the 
fitted parameters of the chiral lagrangian and the inverse 
range parameters aj . The last rows in the tables compare 
our values with those determined in Ref. [4]. We remind 
the reader that the parameter 6o and the baryon mass in 
the chiral limit were not fitted to the data and are given in 
the second and third column of Table [2] only to visualize 
their respective values corresponding to the selected ct^at 
term. The ttN isospin-even scattering length shown 
in the fourth column of Table [2] was not included in our 
fits either but we feel that its presentation is important 
and deserves some comments. 

The goal of the present work was to check the compat- 
ibility of the DEAR kaonic hydrogen data with the low 
energy K~p cross sections and branching ratios. There- 
fore, we have not included in our fits the yl(1405) mass 
spectrum and other processes considered e.g. in Ref. 8 . 
In fact, the low energy constants involved in the fits should 
be also constrained by other observables calculated within 
the framework of ChPT involving the same meson-baryon 
lagrangian. The spectrum of baryon masses and the ttN 
isospin-even scattering length may come to one's mind in 
this respect. The later quantity to order q'^ is given by 

m-. 



Table 3. Inverse range parameters aj (in MeV): 



j+ - 



1 



47r(H-TO^/MAr) 

n2 



f (-2fei)-26F-4feo 



2 

9a 



dD+dF + 2do) - 



(7) 



Since the experimental value of is practically consis- 
tent with zero, = -(0.25 ± 0.49) • m-^ [17], it 
is encouraging to note the mostly negative signs of the d- 
parameters that cancel the positive contributions due to 
the b terms and the q'^ correction represented by the last 
term in Eq. ([7]). As a smaller ct^tv term means a smaller 
absolute value of the negative parameter 6o (and hence a 
smaller positive contribution due to the 6o term in aj^) it 



a^N [MeV] 


CiKN 








Oats 


20 


610 


209 


570 


1100 


530 


30 


647 


262 


535 


308 


21 


40 


653 


320 


618 


281 


89 


50 


594 


370 


610 


342 


124 


27 [5 


760 


300 


450 







is not surprising that the computed -kN scattering length 
is becoming negative for too low ct^at terms. Anyway, it 
is interesting that our fits aimed at the KN interactions 
allow for so good reproduction of the t:N quantity. Specif- 
ically, the parameter set obtained in the fit for a-^N = 30 
MeV gives in a nice agreement with experiment while 
the /N is only slightly inferior to our best fit. Many 
other authors (e.g. |4J or 8J) include the value directly 
in their fits. The d couplings (of the second order chiral 
lagrangian) contribute to the contact meson-baryon inter- 
actions in the second order of meson momenta. Although 
our fits confirm their mostly negative signs it is difficult 
to come to any conclusions concerning their values. The 
fact that even the sign of c?2 is not well determined in our 
analysis speaks for itself. 

We have also tried to perform fits with the b param- 
eters taken from the analysis of the baryon mass spec- 
trum [18] and with only the current algebra (Weinberg- 
Tomozawa) term contributing to the Cij coefficients (the 
approach adopted in Ref. [5J ) . Unfortunately, we were not 
able to achieve satisfactory results in those cases. Thus, 
we conclude that the low energy constants derived in the 
analysis of baryon masses are not suitable in the sector of 
meson-baryon interactions and that the inclusion of the 
q^ terms is necessary for a good description of the KN 
data. The later point is in agreement with the analysis of 
Ref. [7]. 

The inverse range parameters given in Table [3] are in 
line with our expectations. The values corresponding to 
the open channels KN , ttA and ttU seem to be well de- 
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termined and show only a moderate dependence on the 
adopted value of the (Jt^n term. In general, the ranges ob- 
tained for the open channels correspond to the t-channel 
exchanges that are believed to dominate the interactions. 
On the other hand the range of interactions in the closed 
channels is not well defined in the fits and the fitted values 
o^riA/sj and axs exhibit relatively large statistical errors. 
This feature also justifies our use of only one range pa- 
rameter for both rj channels. 

In Figure 1 we present the low energy K~p initiated 
cross sections calculated using our best fit with a^^N — 20 
MeV. The results obtained for the other adopted values 
of a^rN are quite similar, therefor we decided to not in- 
clude them in the figure. Though we declined from using 
all experimental data in our fits and took only the data 
points available for the selected kaon laboratory momenta 
Plab — 110 MeV and plab = 200 MeV, the description 
of the data is quite good. Specifically, we do not observe 
the lowering of the calculated cross sections in the elas- 
tic K~p channel reported by Borasoy et al. |7j for their 
fits including the kaonic hydrogen characteristics. Though 
our K^p cross sections are also slightly below the exper- 
imental data the difference is not significant. In addition, 
the inclusion of electromagnetic corrections discussed in 
Ref. [7] should partly improve the description for the low- 
est kaon momenta. 

Finally, let us turn our attention to the calculated char- 
acteristics of the Is level in kaonic hydrogen. The strong 
interaction energy shift of the Is level in kaonic hydrogen 
is reproduced well but we were not able to get a satis- 
factory fit of the Is level energy width as our results are 
significantly larger than the experimental value. This re- 
sult is in line with the conclusions reached by Borasoy, 
Meissner and Nissler [11] on the basis of their comprehen- 
sive analysis of the K~p scattering length from scattering 
experiments. However, when considering the interval of 
three standard deviations and also the older KEK results 
[10] (which give less precise but larger width) we cannot 
conclude that kaonic hydrogen measurements contradict 
the other low energy KN data. 

In Table S] we compare our results (for a^^N = 20 
and 50 MeV) for the Is level characteristics in kaonic hy- 
drogen with the approximate values determined from the 
K~p scattering lengths aj^-p. The later quantity is ob- 
tained from the multiple channel calculation that uses the 
same parametrization of the strong interaction potential, 
Eq. ID). The Is level complex energies are shown for: the 
standard Deser-Trueman formula (DT) \2:, the modified 
Desert- Trueman formula (MDT) |3j (see also Ref. [7] for 
the relations used to obtain the DT and MDT values) and 
our "exact" solution of the bound state problem. The re- 
sults obtained for ct^a? = 30 and 40 MeV are quite similar 
(with almost identical K~p scattering lengths), so we did 
not include them in the table in order to keep the presen- 
tation more transparent. Since the numerical precision of 
determining the bound state energy by our method is bet- 
ter than 0.1 eV, the discrepancy between the "MDT" and 
the "exact" values should be attributed to higher order 
corrections not considered in the derivation of MDT. In 



Table 4. Precision of the Deser-Trueman formula. The com- 
plex energies AEn — (i/2)/' are given for the approximate DT 
and MDT formulas and compared with our computed "exact" 
values. 







AEn 


- (i/2)r [eV] 




DT 


207 


- (i/2) 832 


-0.50 -l-i 1.01 


MDT 


251 


- (i/2) 714 




exact 


232 


- (i/2) 725 




DT 


247 


- (i/2) 830 


-0.60 + i 1.01 


MDT 


285 


- (i/2) 689 




exact 


266 


- (i/2) 708 



fact, the correction due to Coulomb interaction is taken 
only in its leading order in the MDT formula. The inclu- 
sion of more terms of the relevant geometric series would 
bring the MDT value into a better agreement with our 
exact solution [19] . 



4 Conclusions 

We have computed the characteristics of kaonic hydro- 
gen exactly and compared the results (the Is level energy 
shift and width) with the values determined from the K~p 
scattering length by means of using the standard Deser- 
Trueman formula and its modified version which includes 
the corrections due to electromagnetic effects. It looks 
that the approximate DT formula gives the Is energy level 
strong interaction shift and width about 10% and 15% off 
the exactly computed values, respectively. Although the 
modified DT formula does much better job on account of 
the width the energy level shift remains about 10% off 
the exact value that lies approximately in the middle be- 
tween the DT and MDT values. In view of the current 
level of the experimental precision the use of the modified 
DT formula is sufficient. Nevertheless, the situation may 
change after the coming SIDDHARTA experiment being 
prepared in Frascati. 

An effective chirally motivated separable potential was 
used in simultaneous fits of the low energy K~p cross sec- 
tions, the threshold branching ratios and the character- 
istics of kaonic hydrogen. The fits are quite satisfactory 
except the Is level energy width being much larger than 
the experimental value. In view of the fact that the experi- 
mental precision of the kaonic hydrogen data is still rather 
low one cannot say that the data contradict the chirally 
motivated model used to describe the low energy meson- 
baryon interactions. However, as the opposite statement 
cannot be made either we should wait for the new exper- 
iment to clarify the situation. 
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Fig. 1. Total cross sections for K p scattering and reactions to the meson-baryon channels open at low kaon laboratory 
momenta plab- The experimental data are the same as those compiled in Fig. 1 of Ref. [4]. 
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